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Алгоритм для решения нелинейных 


краевых задач по т-методу Ланцоша 
в системах компьютерной алгебры 


Построен алгебраический алгоритм для преобразования многоточечной линейной краевой задачи для 
дифференциального уравнения порядка К с линейной частью — линейный дифференциальный оператор 
многочленными коэффициентами порядка А и нелинейной частью — функция Ху, у’...., УГ) в 
алгебраический многочлен порядка и е №. Этот многочлен — аппроксимация решения у(х), х е [а,6] 
исходной краевой задачи. Эта аппроксимация оптимальна в пространстве Сы: 


Введение 


Актуальность темы исследования. Многоточечные краевые задачи для обык- 
новенных дифференциальных уравнений являются классическим аппаратом матема- 
тического моделирования [1]. Маре, Мафетайса, Масаа, МаНаБ и другие системы 
компьютерной алгебры (СКА) — системы символьных вычислений на компьютерах — 
стали естественной средой математического моделирования. 

СКА для обыкновенных дифференциальных уравнений символьно вычисляют: 

— решение — композицию специальных математических функций (существует редко), 
— частную сумму ряда Тейлора решения задачи Коши порядка п. Обычно п < 10. 
Этот алгебраический многочлен, как правило, не удовлетворяет основной критерий 
эффективности математических моделей — точность. 

Возможности классического математического обеспечения ЭВМ. Компьютерные 
системы типа Ма Наб вычисляют сеточную аппроксимацию решения краевых задач для 
обыкновенных дифференциальных уравнений по численным методам. Это методы, как 
правило, с насыщением. Производные сеточной функции, как правило, являются не 
оптимальным аппаратом аппроксимации производных исходной функции. 

Цель работы. Построить алгебраический алгоритм для вычисления в СКА ап- 
проксимации на отрезке [а,Б] решения многоточечных линейных краевых задач для 
обыкновенных дифференциальных уравнений порядка А следующего вида 


БЫ =ЛУ, у... У,  Шу]=А +... +Су, (0 

Сота») = { РДУ] в -а:+ б=0}:-и, БУ = АУ +... + Су. (2) 

Коэффициенты 4, ....С линейного оператора О[у] (1) — алгебраические многочлены. 

Точки задания краевых условий (2) принадлежат отрезку [а,6] — 4 Е [а,6], 1=1,..., К 
и являются регулярными не особыми точками оператора О[у]| - А(4;) = 0,1=1, ..., К. 
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По этому алгоритму в СКА вычисляют алгебраический многочлен у„ порядка п Е № 
У=с+ах+... Неих" = авотйт_1((Т), (2), [а,6], п). (3) 


Многочлен у, (3) аппроксимирует решение у = 50е( (1), (2) ), у = у(х), х Е [а,Б] 
краевой задачи (1), (2). Эта аппроксимация оптимальна для преобразований в СКА 
решения у задачи (1), (2) и его производных у’. ..., У® на отрезке [4,6] — коэффициент 
оптимальности алгоритма на краевой задаче (1), (2) в пространстве С’ми,ь] ограничен 
С(авоптйт_1, (1), (2), Сы) = [у- у» [х/ МЕ 0, се» | У) - с+еах+...+ сих”) |х <. 

Актуальность задачи. Уравнения вида (1) применяются в математических мо- 
делях более часто, чем дифференциальные уравнения других типов. Если процесс 
моделируют дифференциальным уравнением порядка (А, то, часто, на отрезке [а,Ь] 
исследуют решение у = у(х) этого уравнения и его производные у’ ..., у — преобра- 
зуют эти функции на отрезке [а,Б]. 


1. Композиция итерации, т-метода и интерполяции 


Алгоритм 1. 

Вход: Уравнение (1), условия (2), отрезок [а,6 |, параметр и. 

Выход: Алгебраический многочлен (3) (порядка п ). 

Преобразования: 

1. Вычислить начальное приближение к решению краевой задачи (1), (2) — мно- 
гочлен порядка А — 1, удовлетворяющий краевые условия (2) 


Уно = Ут = 50[е( Сопа( ут Е Ри!) ). (4) 
2. Для5= 1,2, ... ВЫПОЛНИТЬ: 
2.1. Вычислить производные у’ ;-1, ..., а многочлена ув, ;-1. 
2.2. Преобразовать функцией / (1) многочлен ул,;.1 и его производные в функцию 
О но (5) 
2.3. Вычислить алгебраический многочлен 
Е; = ОК] = 94 (а + (Б-а) (1- сот /п ))/2), 1=0,..., п}] - (6) 


интерполяцию функции / (5) в узлах Чебышева на отрезке [а,в] [2, с. 94]. 
2.4. Вычислить линейное дифференциальное уравнение с многочленными 
коэффициентами (ЛДУМК) 


РЫ = Ро. (7) 


2.5. Решить по алгоритму [3] х-метода Ланцоша [1] с параметром п краевую 
задачу (2) для ЛДУМК (7) на отрезке [а,6] и вычислить аппроксимацию точного ре- 
шения краевой задачи (7), (2) — у(5,х) = 50[»е( (7), (2) ) — алгебраический многочлен 
порядка и 

у»; = авогийт_[3]( БЫ] = Е, Сопа(»), [а,6], п). (8) 


3. Вычислить искомую аппроксимацию (3) решения у краевой задачи (1), (2) — 
предел последовательности многочленов (8) при 5 —> © 


у = а[еогийт_1 ( (1), (2), [а,6], п) = Ша; У. 
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2. Оптимальность алгоритма 1 на модельной задаче 
Модельная краевая задача метода квазилинеаризации [4] 


У" = ехро), 0 )=0, (1) =0. (9) 
Эта задача является одномерным аналогом уравнения гидродинамики у", + у", = ехр(») 
и имеет единственное решение — функцию аналитическую в окрестности отрезка [0,1] 
комплексной плоскости 
у=- (2) + 2 ш(с зес(с (х- 1/2)/2), с = зо»е(за (2) = с зес(с/4)). (10) 
Свойства функции у (10). Отличные от нуля коэффициенты Фурье — Чебышева 
функций у, у" (10) на отрезке [0,1] только четные и принимают следующие значения 
{ ак», [0,1]) 1-0 = {$ -0.057, 0.057, 0.00027 , 2.1 10°, 1.8 10%, 
1.66 10°, 1.6 10^?, 1.6 107“, 1.6 107°, 2.1 1078}, 
{ аку", [0,1]) };-0°= $0.94, 0.053, 0.001, 1.6 102, 2.4 107, 
3.5 10°, 4.71071, 6.3 10°, 8.2 10°, 8.8 1077, -5.5 1078}. (11) 
С ростом параметра и = 2 Тони регулярно убывают к нулю 


[42 1+2(У, [0,1] / [42 Ку, [0,1] = у, 2ь ую = 1, 1,2 = 0.005, 1,4 = 0.007, эк» УТ: 0.01, 
[42:+2(у’,, [0,1] / [42 КУ”, [0,1] = 2,2 20 = 0.06, 22 = 0.2, 24 = 0.016, 4:4 /2,16 = 0.01. 
Следовательно, для величины наилучшего приближения функции у (10) алгебраичес- 
кими многочленами порядка и в пространстве Ср1] и Сол] справедливы тождества 


Ш 0, дел || У) — (©... + их”) [сю = (1 + В») [ 42 21 +20, [0,1] |, 
Ш 0, дея || 709) — (©+... Ни х”) [с:2 вп = (1 + 9) | а2 из”, [0,1]) |) — (12) 
гдеи> 1 - | В„ | < В= 0.01, и>2 — | 7, | <у= 0.01. 
Вычислительный эксперимент с алгоритмом 1. Норма погрешности решения 
краевой задачи (9) по алгоритму 1 в пространстве Су] и Со] удовлетворяет тож- 
дества 


{|у- у: [ео г = { 0.0043, 2.5 10°, 1.2 107, 1.1 10°, 6.7 107! }, 
{у-у2; [сю иг = { 0.056, 0.001, 1.7 10°, 2.5 107, 3.6 107}. 


Из тождеств (11) и (12) можно сделать следующее заключение. 
Вывод 1. Главной частью коэффициента оптимальности решения по алгорит- 
- 2 
му 1 краевой задачи (9) в пространствах Сю,ци С`юл] являются отношения 


{Пу У» Паола / | 42 изу +20% [0,1]) |= С›изь $ Со" = { 16, 12,7, 6 }, 
[у — Ун||с^2_[0, 1/2 [п/2] ("0,11 =1+ [9.2% [п/2] т — {0.06, 0.000006, 0.05, 0.04}.(13) 


Согласно определению, коэффициент оптимальности метода вычисления алгебраи- 
ческого многочлена, аппроксимирующего решение задачи (1), (2), не меньше единицы. 
Поэтому из этих тождеств и тождеств (12) можно сделать следующее заключение. 

Вывод 2. Главная часть коэффициента оптимальности алгоритма 1 на зада- 
че (9): 

— в пространстве Су] ограничена и справедливо тождество 


С; авогийт_1, (9), Сон) = (1+ В») С>ш 


— в пространстве Со асимптотически тождественна минимальному значению коэф- 
фициента оптимальности и справедливо тождество 


С‚( авогийт_1, (9), Сюц) = (1+ у,) (1+ @2 р] ). 
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3. Алгоритм 1[ и метод квазилинеаризации 


Метод квазилинеаризации [4]. По этому методу решают двухточечные краевые 
задачи для нелинейных дифференциальных уравнений второго порядка, разрешен- 
ных относительно старшей производной (1а45^) с точками задания краевых условий а 
и В (типа задачи (9)) на отрезке [а,6], и вычисляют сеточную функцию с 4 + 1 равно- 
отстоящими узлами сетки 


х=а м=а-й,....жм=Ь, й=(Б-а)/4, 
ул = у 0" = $ ул (к, ак), -07 = тефо4[4]( а5к, [а,6], а) = УС". 


Композиция метода квазилинеаризации и интерполяции алгебраическим мно- 
гочленом сеточной функции у» (вычисленной по методу квазилинеаризации) является 
методом вычисления алгебраического многочлена порядка 4 


Ра у» | = Р[тейо4[4]( 1а5к, [а,6], 4)], й = (Б-а)/4. (14) 


Этот многочлен аппроксимирует решение у исходной задачи {5 на отрезке [а,6]. Норма 
погрешности метода (14) в пространстве С|ь] не меньше соответствующей сеточной 
нормы. Сеточная норма погрешности метода (14) тождественна сеточной норме по- 
грешности метода квазилинеаризации 


ах, с [вы [У-— Ра[У+] | > Шах; -0,....а| -— Ра] к =х+| = тах; 0, ... а| У(х)) — ух |. 


Следовательно, коэффициент оптимальности метода (14) в пространстве Х= Сщиы не 
меньше 


ССРаун 195%, Сы ) > тах: =0,...а| Уса) — УК) | / Ш о,..е а || У) — (с... + см) [са 


Сравнение алгоритма 1 с методом Р.[у»]| (14). На краевой задаче (9), согласно [4, 
с. 42], метод квазилинеаризации с параметром 4 = 10 имеет погрешность 3 10’. Поэтому 
из тождеств (11) — (13) можно сделать следующее заключение. 

Вывод 3. Главная часть коэффициента оптимальности в пространстве Суо,11 
метода алгебраической интерполяции сеточного решения метода квазилинеариза- 
ции (14) на краевой задаче (9) более чем в 100 000 раз больше 1 + © шир] (13) — главной 
части коэффициента оптимальности алгоритма [ на этой задаче 


Сто( Раоуь = мо] › ва8К (9), Сюи)>3 107 1.6 107? > 10°. 


4. Алгоритм 1 и методы с насыщением 


Метод квазилинеаризации и другие сеточные методы решения функциональ- 
ных уравнений, как правило, являются методами с насыщением. Скорость убывания 
с ростом числа узлов погрешностей решения по методу с насыщением краевой зада- 
чи с решением — аналитической функцией полиномиальная 


тах;-0, .. | (д — ух = О(а?), й=(Ь-а)/ 4, р = р(тейоа) = Сопз4. 


Скорость убывания (с возрастанием параметра п) величины наилучшего приближе- 
ния функции у (10) алгебраическими многочленами порядка п в пространствах Суо,11 
и Со] оценивают тождества (11). Поэтому коэффициент оптимальности метода (14) 
и других методов интерполяции сеточного решения у„ многочленами возрастает с рос- 
том числа узлов сетки 4 и скорость роста 


ССР, 1А5К, Сю) = 0(2'), ай 
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Коэффициент оптимальности метода ряда Тейлора имеет такую же скорость роста. 
Из этих тождеств и вывода | можно сделать следующее заключение. 

Вывод 4. На краевой задаче (9) алгоритм 1 предпочтительнее метода (14) и 
метода ряда Тейлора, согласно критерию — минимум коэффициента оптимальности. 


5. Эквивалентный алгоритм метода Галеркина 


Пусть для условий (2) существуют многочлен ух_ 1 (4) и линейный интегральный 
оператор 


КИ] = золе( У = { БУ а =0 1-1"), (15) 
где О} у | - линейные дифференциальные операторы условий (2). Тогда результат под- 
становки у = КГ и ] + ук _1 в уравнении (1) является интегральным уравнением отно- 
сительно функции и 

В КТи +1] =Д КГ и] + ть СКТи] + ук) ..› (КГи Ни). (6) 
Для решения и уравнения (16) и решения у краевой задачи (1), (2) справедливо тож- 
дество 

У=КГ и] + ук-1. 

Метод 1. 

1. Преобразовать краевую задачу (1), (2) в интегральное уравнение (16). 

2. Вычислить аппроксимацию уравнения (16) по методу Галеркина с оператором 
проектирования 5»: [2(а,Б;р) —> Н‚а,Б] (вычисляет частную сумму порядка р = п - ^, 
где А — порядок уравнения (1), ряда Фурье — Чебышева функции на отрезке [а,6]) и 
возмущением алгебраическим интерполированием правой части уравнения (16) по уз- 
лам Чебышева на отрезке [а,6] — („[./] (6) 


З»ГРЕ Кир 4-1 1 = ЗС ОСА Кир Нет... Кир к-т)" (и © НЫа,Б.(1Т) 
Неизвестным уравнения (17) является алгебраический многочлен и, порядка 
р =п- вида — частная сумма ряда Чебышева на отрезке [а,6] — и, е Н]а,Ь]. 
3. Решить уравнение (17) по методу простой итерации с начальным приближе- 
нием иро = 0. Следующие приближения являются решением уравнения 


Ир = 50№е( 5 Б[ К и | +у-11-6]=0 (ше НЫа,Б] ) ), (18) 
а 


По этому методу вычисляют алгебраический многочлен и, = Шт; Ир». 
4. Преобразовать многочлен и, в искомую аппроксимацию решения у задачи (1), (2) 


У» = КГир ] + Ук-1. 

Замечание 1. Уравнение (18) является аппроксимацией по методу Галеркина ли- 

нейного интегрального уравнения с многочленными коэффициентами (ЛИУМК) 
ВК и] + ук, =А.. (19) 

Уравнение (19) является результатом подстановки у = К|[ и ] + ук! (введена в уравне- 
нии (16)) в ЛДУМК (7). Следовательно, метод 1 можно реализовать в виде следую- 
щего алгоритма. 

Алгоритм 2. Алгоритм 1 со следующей модификацией п. | ип. 2.5. 

1. Вычислить начальное приближение к решению уравнения (17) и» = 0, опе- 
ратор К (15) и многочлен (4) у = У (КТО ]=0). 

2.5. Вычислить: 
— ЛИУМК (19) — преобразовать краевую задачу для ЛДУМК (7), (2), 
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— решение и,,; (18) уравнения (19) по методу Галеркина, 
— аппроксимацию решения задачи (7), (2) — преобразовать многочлен ир,; (18) 
Уи, = К] Ир, Пе ыяа (20) 

Теорема 1. Пусть: 
— для условий (2) существуют многочлен у. (4) и оператор К] и ] (15), 
— оператор К[] и ] (15) преобразует моном х, 1 Е М в полином 

Кх] = хе + ‚аеь В; = 0, 

— параметр алгоритма 1п> 2, К= от еди( Бу] ) - порядок уравнения (1). 

Тогда алгоритм 1 эквивалентен алгоритму 2. 

Доказательство. Если выполнены условия теоремы, то, согласно [3], многочлен 
У„з (8) тождественен многочлену у»; (20). 

Остальные преобразования алгоритмов 1 и 2 тождественны. 


6. Сходимость алгоритма 1 


Из оценок оптимальности операторов 5, и (И, как аппарата аппроксимации функ- 
ций на отрезке [а,Б] [2, с. 77-95] 


Пу [У 1 25) = о, ет [| У(Х) — (0+... +) [12.5 
|5» [ыы = 1, || 5 аьы = (4/2) № п) + 0), 0) $3, п>0, 
ТО, [авы = (2/5) Сп) + О(0), ОО <1, п>0 
можно сделать следующее заключение. 

Вывод 5. Результаты исследования метода Галеркина [5] решения операторных 
‘уравнений в пространстве Гильберта [2(а,Б; р) на интегральных уравнениях вида (16) 
достаточно широкого класса доказывают: 

— существование многочлена 


ир= зоие( (17)), рт, т=т( БК и], Ажуу, УК мини ), 
— сходимость последовательности Иш, ..., Ир, Иры, ... — и = 50№е( (16)), 
— ограниченность коэффициента оптимальности метода 1 (и алгоритма 2). 
Поэтому из теоремы 1 для краевой задачи для обыкновенных дифференциаль- 
ных уравнений вида (1), (2) достаточно широкого класса следует: 
— существование многочлена у» (3), п > т-+к, 
— сходимость последовательности ут, ...› Ув Уныь ... =» У = 5О№е( (1), (2) ). 


7. Апостериорные оценки погрешности 


Теоретические оценки. Для алгоритма [3] известны точные и конструктивные 
апостериорные оценки нормы в пространстве Ста] и СЕ 5] погрешности решения 
(5, х) краевой задачи для ЛДУМК (7), (2) — || У(5,х) — у„(х) ||. Если последовательность 
этих оценок (5 = 1, 2, ...) сходится, то ее предел естественно принять за оценку по- 
грешности алгоритма 1 на краевой задаче (1), (2). 

Пример 1. (Иллюстрация эффективности теоретических оценок). Мы решили по 
алгоритму | краевую задачу (9) и тождественную задачу 


У-у=ехро)-у, У(0)=0, УП =0. (21) 
Начальное приближение у.о = у! = 0. По алгоритму 1 для задачи (21) решают сле- 
дующие краевые задачи для ЛДУМК 


уу = Ц) ехрО»» 5-1) — Ул 51|, У(0)=0, УГ) =0, 5=1,2, .... (22) 
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Оценка [3] нормы в пространстве Х погрешности решения по алгоритму [3] краевой 
задачи (22) на отрезке [0,1] с параметром и имеет вид 


|| И (и/21+2(%) ||х | (в $), 5=1,2, ..., (23) 


где т(и,5) = <(2 [п/2], 5) — отличный от нуля коэффициент дополнительного многочле- 
на решения задачи (22) по алгоритму [3], 


| И > ||с2 юи=1, 1=0,1,.... 
Коэффициенты { *(6,1), (6,2), ...} оценки (23) принимают следующие значения 
{ 610’, 1.7 10?, 1.65326 10°, 1.65277 107,...,1.65277 10°,...}. 


Последовательность т(и, 1), т(и, 2), ... (пЕ №) коэффициентов оценки (31) имеет пре- 
дел (и) = т(2 [п/2]). Этот предел для задач (9) и (21) принимает следующие значения 
{ (2), т(4),...2(10) } = { 0.056, 0.001, 1.7102, 2.5 10”, 3.5 10? }. 
Поэтому из оценок [3] и тождеств (11) - (13) можно сделать следующее заключение. 

Вывод 6. /Гредел оценок [3] погрешностей алгоритма [3] на краевой задаче для 


ЛДУМК (7), (2) при 5 —> © является конструктивной оценкой погрешности алгорит- 
ма 1 на краевой задаче (1), (2). На краевой задаче (9) и (21) эта оценка точная. 


8. Априорные оценки погрешности 


Теоретические оценки. Для решения краевой задачи для ЛДУМК (7), (2) по ал- 
горитму [3] известны [3] точные и конструктивные априорные оценки нормы погрешнос- 
ти и коэффициента оптимальности в пространстве С* [а, 5. Этот коэффициент зависит 
от оператора БГ у ] (1) и операторов ДР; (2). Если линейная часть (по у) функции Г 
уравнения (1) мала, то коэффициент оптимальности решения краевой задачи для 
ЛДУМК (7), (2) по алгоритму [3] достаточно точно аппроксимирует коэффициент оп- 
тимальности решения краевой задачи (1), (2) по алгоритму 1. 

Пример 2. (Иллюстрация эффективности теоретических оценок). Оценка [3] 
коэффициента оптимальности в пространстве Сон] решения по алгоритму [3] задачи 
(22) на отрезке [0,1] имеет вид 

Ск(авогийт_[3}, (22), С’ю, п ) = || И" К® сю, п/а2 К И" Кх), [0,1] ), (24) 
где 1= [и/2]. Мы вычислили составляющие правой части оценки (24) — 
|” Кх) [<ю.и =1, 1/2 КИ”> Кх), [0,1]) =1+ 62ь 
{ 02,1=0, ...,6 } = { 0,06, 0,04, 0,008, 0,0036, 0,0022, 0,0013, 0,00087 }. 
Из этих тождеств и тождеств (13) можно сделать следующее заключение. 
Вывод 7. В пространстве С? Го, 11 Оценка (24) коэффициента оптимальности алго- 


ритма [3] на краевой задаче для ЛДУМК (22) является эффективной оценкой коэф- 
фициента оптимальности (13) алгоритма 1 на краевых задачах (9) и (21). 


9. Вычисление оценок погрешности 


В работе [3] построены следующие алгоритмы для вычисления оценок погреш- 
ности решения краевых задач (2) для ЛДУМК (7) по алгоритму [3]: 
— алгоритм 2 оценки нормы погрешности в пространстве Си в, ..., 
— алгоритм 3 оценки коэффициента оптимальности в пространстве Са, 5]. 
Для задачи (22) по этим алгоритмам с параметром 4 вычисляют многочлен 


И’. а(х) = И’ (х), || 6) — И) [с то, и < (и, 4) 
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и один отличный от нуля коэффициент дополнительного многочлена х(и, 4). Этот ко- 
эффициент убывает к нулю (с ростом параметра 4 ) 


х(п, 4) |= о(и1`"), и< 0.1. 
Поведение коэффициента *(и, 4) хорошо иллюстрируют его значения 
{2 10) -о0’= {2 1076,3 107°, 8 1073,4 107°, 3 107, 5 10“, 1}. 
Из этих тождеств можно сделать следующее заключение. 


Вывод 8. В случае задачи (22), отрезка [0, 1] и параметра п: 
— по алгоритму 2 [3] с параметром 4 вычисляют аппроксимацию 


|| #7 изн, а(%) [|с^2_по, п [| <(м, 5), 5=1,2,..., 
оценки (24) погрешности алгоритма 1 и погрешность этой аппроксимации не больше 
ИР [2+ 26%) [2 101 — [| 2 ин, 40%) [| с^2 до | | (и, 5) | < 
| ч(и, 5) | (2 [и/2] + 2, 9) [= [тб $) [об ОР) = о), и< 0, 
— по алгоритму 3 [3] с параметром 4 вычисляют аппроксимацию 
| И [и/2 +2, (%) ||сто,11 а? [2 +2(Й”5 м/о, а(%),[0,1]), 

оценки (24) коэффициента оптимальности алгоритма [ и погрешность этой аппрок- 
симации имеет порядок 

О(|т(2 [и/2]-+2, 9) |) = Ой?" ), и<0.1. 


10. Программирование алгоритма 1 в АР 


Структура данных на входе. 

1. Уравнение (1) определяют его правая и левая части — функция Ху, у’,..., У 
и линейный дифференциальный оператор с многочленными коэффициентами Б]| у |. 

2. Оператор В[ у] имеет вид А * О1Е( у, К) +... + С * у + С, 
где у -атом, многочлены А,..., С, С являются термами, имеют вид, естествен- 
ный для математики и аргумент -— атом х. 

3. Функция Хх, у, у, ...,Э) имеет вид, обычный для СКА, и аргументы — атомы 

к УЗ Уря УБЕЙ 

4. Условия (2) имеют вид (...,(х = 4 1, А*аЩу,К_1)+ ... +С*у-+Сп = 0)....). 

5. Отрезок аппроксимации [а, Б] определяет список (а, Ъ). 

6. Параметр и алгоритма является целым числом. 

Структура данных на выходе. Многочлен у, (3) имеет коэффициенты -— числа 
и вид, естественный для математики, а +... + Ё * х^п. 

АРГАМ-процедура. Алгебраическая спецификация п. 2 алгоритма 1. 


уг а = ахум; а Е 
у2_п := Я х(у1 п); И Е. 2 
И ве У р УЛ РО УРА У ДЬ 
вы = сне тобетро(Еуче, пабе а п ав 
ЪООМК := (Пу (-1) * Ез = ЕЯ 


АРГАМ-процедура - реализация алгоритма [3]. 


Структура выхода операторов АРГАМ-процедуры. 
Процедура вычисляет начальное приближение уд (4) многочлен с числовыми 
коэффициентами вида, естественного для математики. 
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На каждой итерации (5) процедура вычисляет: 
— производные многочлена у„;.1 вида, естественного для математики, 
— функцию =ХКх, уу...) |= 5-1}; естественного вида, 
— многочлен А, = („| ] естественного вида, 
— ЛДУМК (7) вида [3] А*АТЕ(у, К) +... + С*у+ -1 * Ез=О0, 
— следующее приближение у»; к решению исходной краевой задачи (1), (2) — много- 
член с числовыми коэффициентами вида, естественного для математики. 


11. Исследование АРГАК-процедуры 


Теорема 2. Сложность выполнения одной итерации алгоритма 1 АРГАМ-про- 

цедурой с ростом параметра п возрастает как полином 
(т-+1) О( сапрИ, т) + О( п’), т= шах { 4е2( Бу, еЕР,] ), п} + =п+ О(\), 

где О( сапр!, т ) — сложность преобразования оператором сапр1Е многочлена по- 
рядка т к каноническому виду — сумме мономов вида с (1) *х^) $5, 1,1=0, ..., т. 

Доказательство. Процедура вычисления одной итерации решения линейная. По- 
этому вычислительная сложность этой процедуры тождественна сумме вычислительной 
сложности операторов этой процедуры. Вычислительная сложность АРГАМ-процеду- 
ры, реализующей алгоритм [3], оценена в работе [3]. Остальные операторы процедуры 
имеют по параметру и полиномиальную сложность О(и”), $ < 3. 


12. Вычислительный эксперимент с процедурой 


Описание краевой задачи (21) на языке АРГАМ. 


ргосез$[1] := ( ру = ЕС 2 ему; 
Еу := ехр(у) + -1 * у; 
бопа = (= брус =о)): 
1п6егуа1 := (0,1); 2); 
Результаты решения задачи (21) процедурой (и = 2). 
Ру: 3, ВОВ уе. СУ Ме Ле 
Е. :=.спев. таветрот(Еу м „ пивеетат в = 1; 
ПОМИ: вес № РВ =0) = 
(ЧЕ (у, 2) + -1 * у+-1=0); 
Зы =, БОТОВ) <=10. 410588: 5” 2+ =0:410588 5.3 
К 5 := СнеБ. тпбетро1 (Еу п ; лабетуат , п} = 


—0.0266273 *х ^ 2 + 0.0266273 *х+1; 

ЪРОМК := (Ру + -1* Евз=О) = (Е (у, 2) + -1 ху+ 
00266273: * х^2+ 0.026623 хо; 

у о верисовЕ) = оао. х Пр ох 

У ее СОвЕ) = беда О 5: Е. Е а 65 2} 


13. Программирование алгоритма 1 в СКА 


Построенная выше АРГАМ-процедура доказывает — алгоритм 1 имеет алгебраи- 
ческие преобразования и его можно реализовать в других СКА. Результаты исследо- 
вания этой процедуры доказывают эффективность такой реализации. 


«Штучний 1нтелект» 472009 127 


Денисенко П.Н. 


ЗД 


14. Решение практической задачи — задачи Прагера 


Задача о деформации упругой струны под действием поперечной нагрузки. 
Конечные деформации упругой струны под действием поперечной нагрузки хорошо 
моделирует краевая задача 

у"=1+а” (у, У0)=0, У) =0. (25) 
Эту задачу предложил В. Прагер для оценки эффективности метода квазилинеариза- 
ции [4, с. 43]. Задача (25) имеет единственное решение — функцию, аналитическую в 
окрестности отрезка _ П комплексной плоскости 

= Иа” ш( со$( а (х-1/2) ) / соз( а/2 ) ). (26) 

Свойства и у (26). Отличные от нуля коэффициенты Фурье — Чебышева 
функции у (26) и функции у" на отрезке [0,1] только четные. В случае параметра а = 0,7 
они принимают следующие значения 

{ а2 К», о 0) - О: 064, -0.064 , -0.00034 , -2.8 107 И 
7101, -2.9 1072, -3.2 10°, -4.6 10-16, 6.2 1018}, 
2 ", [0,1] };- в - о т -0.067 , -0.0014 , 2.5 10°, -4107, 

-6.3 103, -9.5107", -14 1072, 2107“, -2.6 1076, 4.41077}. (27) 
С возрастанием параметра п= -2 Гони регулярно бывают к нулю, аналогично функ- 
ции (10). Следовательно, для величины наилучшего приближения функции у (26) с 
параметром а = 0.7 алгебраическими многочленами порядка и в пространстве Сю 1] и 
в пространстве Сю, справедливы аналоги тождеств (12). 

Вычислительный эксперимент с алгоритмом 1. Норма в пространстве Со] 
иС [0,1] погрешности решения по алгоритму 1 задачи (25) с параметром а = 0,7 удов- 
летворяет тождества 

{Пу-у2; Сюй г = = {0.006 , 3.7 102, 1.8107, 1.8 10-9, 4.8101}, 
у-У2: сю Нет = { 0.07, 0. 0014, 2.5 10°, 4.1107, 10% }. (28) 
Из тождеств (27) и (28) можно сделать следующее заключение. 

Вывод 9. Главной частью коэффициента оптимальности решения по алгорит- 

му 1 краевой задачи (25) в пространствах Сул] и С’ю,ц являются отношения 
ПУ У» [©1011 / | 42 1 +2(у › [0,1]) |= С? Сы] =1° = { 18, 13, 6.(6), 6.(6), 6.2 }, 
|’ Ус юа2 [ п/2 10", [0, Ц = =1+ 02 [п/2]> {09 =| = = { 0. 05, 0. 009, 0. 004, 0.002 р (29) 

Согласно определению, и оптимальности метода вычисления алгебраи- 
ческого многочлена, аппроксимирующего решение краевой задачи (1), (2), не меньше 
единицы. Поэтому, из этих тождеств и тождества (12) можно сделать следующее зак- 
лючение. 

Вывод 10. Главная часть коэффициента оптимальности решения по алгорит- 
му 1 краевой задачи (25) с параметром а = 0.7: 

— в пространстве Сул] ограничена и справедливо тождество 

С @еотийт_Т, (25), а= 0.7, Сп ) = (1 ыы В,) С м2» > 1 В, < В = 0.015, 

— в пространстве Сол] асимптотически тождественна минимальному значению 
коэффициента оптимальности и справедливо тождество 
С @еотийт_1, (36), а= 0,7, Со.) а (1 + у») (1 + @2 [2/2], п>2-> Уп <= 0.015. 

Сравнение алгоритма 1 и метода квазилинеаризации на задаче Прагера. 

На краевой задаче (25) с параметром а = 0.7, согласно [4, с. 42], метод квазилинеа- 
ризации имеет погрешность 3 т Поэтому из тождеств (27) - (29) можно сделать сле- 
дующее заключение. 

Вывод 11. На краевой задаче (25) (а = 0.7) главная часть коэффициента опти- 
мальности метода алгебраической интерполяции сеточного решения метода квази- 
линеаризации (14) (4 = 10) более чем в 1 000 000 раз больше 1 + 02 [р] (29) — главной 
части коэффициента оптимальности алгоритма 1 на этой задаче 

Со(Ри[ у» |, ак (25), а= 0.7, Сю.) > 3 10° /2.9 107"? > 10°. 
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Заключение 


Построенный в работе алгоритм 1 имеет только алгебраические преобразования 
и решает весь класс краевых задач для обыкновенных дифференциальных уравнений 
порядка К вида (1), (2) с коэффициентами линейного оператора Б[ у ] (1) — многочле- 
нами и достаточно широким классом функций Хх, у,у'...., У”). По этому алгоритму 
в СКА вычисляют алгебраический многочлен у„ порядка и е М- аппроксимацию реше- 
ния у = у(х),х Е [а, Ь|, исходной краевой задачи (1), (2). Эта аппроксимация оптимальна 
для символьных преобразований функции у(х), х е [4,6] и её производных у’, ..., уУ®- 
в пространстве С\а.ь] ограничен коэффициент оптимальности алгоритма 1. Следова- 
тельно, на краевой задаче (1), (2)) алгоритм 1 реализует идею К. Ланцоша. На основа- 
нии т-метода построить эффективные алгоритмы для вычисления алгебраических 
многочленов — аппроксимации решения функциональных уравнений отдельных типов. 


Дополнение 


1. Алгоритм 1". Модификация алгоритма 1. Алгоритм 1 [3] заменен алгоритмом 1' [3]. 
Это алгоритм т-метода с базисом пространства Ни, ь. 

2. Для алгоритма 1" имеют место аналог теоремы 1 и аналоги применения тео- 
рем [3]. В аналогах пространство /Г2(а,Ь; р) заменено пространством Ни, ь. 

3. Модификация АРГАМ-процедуры. Процедура — реализация алгоритма 1 [3] — 
заменена процедурой — реализацией алгоритма 1' [3]. Она реализует алгоритм 1". 

4. Оценка погрешности алгоритма | в комплексной плоскости следует из тож- 
деств теорем [3]. 

5. Оценка погрешности алгоритма 1" в комплексной плоскости следует из тож- 
деств аналогов теорем [3]. 
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П.М. Денисенко 

Алгоритм для розв’язування нел1и них крайових задач 

за т-методом Ланцоша в системах комп’ютерно! алгебри 

Побудовано алгебрайчний алгоритм для перетворення багатоточково! лийно! крайово! задач! для 
диференщального равняння порядку А з линною частиною — лйний диференщальний оператор з 
коефищентами — многочленами та нелиййною — функщя Ху, у'..., \®?)) на алгебрачний многочлен порядку 
п Е М. Цей многочлен — апроксимащя розв’язку у(х), х е [а,6] орипнально! крайово! задач. Ця апроксимащя 
оптимальна в простор! С\иь] 

Р.М. Реш5епко 

Зо!ушо фе Моштеаг Воипдагу-Уаше РгоМет Ошо Ве Г.апс20$ 

х-Мешод4 ш е Сотршег Асебга Зуетз$ 

Ме сопзитасе4 Фе а1оебгалс а]еоги та Юг гапзЮгите фе поптеаг Боипдагу-уале ргоМет пю Фе а1зебгалс 
ро[упопа1 оё огаег п е М. ТЬ1$ ро]упопиа| 15 фе зо[айоп у(х), х е [а,Б] арргохипайоп Юг Фе ргоМет. 
ТЬ!$ арргохипанов 15 орйта! 1 Фе зрасе С* [,5]- 
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